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Abstract 



' We study the antiferromagnetic Potts model on the Erdos-Renyi random graph. By identifying 

, a suitable interpolation structure and proving an extended variational principle we show that the 

replica symmetric solution is an upper bound for the limiting pressure which can be recovered in 
the framework of Derrida-Ruelle probability cascades. A comparison theorem with a mixed model 
, made of a mean field Potts-antiferromagnet plus a Potts-Sherrington-Kirkpatrick model allows to 

■ show that the replica symmetric solution is exact at high temperatures. 
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^ '. 1 Introduction and main results 

In this paper we give some rigorous results on the antiferromagnetic Potts model on the Erdos-Renyi 
' random graph. 

, It is well known that antiferromagnetic Potts models on graphs are related, at zero temperature, 

' to the graph coloring problem which consists in placing colors on the graph vertices in such a way 

, that two of them connected by an edge have different color. In recent time an algorithmic approach 

was developed to study their colorability [7j based on methods and ideas from disordered system, in 
particular on the replica symmetry breaking scheme introduced within the mean field theory of spin 
glasses ,21 . 

Since the Erdos-Renyi random graph has a locally tree-like structure and large loops, the statistical 
^ • mechanics model with antiferromagnetic interactions has been reported to display some spin glass 

behavior in the physics literature |25) . In particular it has been argued that the one-step replica 
symmetry breaking solution does not get improved by a higher number of steps |28) . 

The rigorous theory of spin glasses has, on the other hand, made important progresses in the last 
decade by means of the Guerra-Toninelli |16j interpolation scheme, their thermodynamic limit control, 
the Guerra-Talagrand |23j theorem for the free energy of the Sherrington-Kirkpatrick model and with 
the general scheme introduced by Aizenman-Sims-Starr [1] which includes an extended variational 
principle to obtain the exact solution for a large class of models. 

This paper is a first attempt to derive within mathematical rigor some properties of the antiferro- 
magnetic Potts model on the Erdos-Renyi random graph. A full treatment of the ferromagnetic Ising 
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case has been given in for focally tree- like random graphs and extended in [TT]. The techniques 
used there are heavily based on the use of ferromagnetic GrifSths-Kelly-Sherman and Grifhths- Hurst- 
Sherman inequalities and do not apply to our case. The extension to Potts random variables is an 
open problem. 

For the antiferromagnetic model we introduce here an interpolation scheme and prove its mono- 
tonicity. It had been previously derived [5] a new discrete version of the Guerra-Toninelli interpolation 
argument to prove existence of the thermodynamic limit for this model for all T > 0, as well as directly 
for T = . This gave a positive solution of a conjecture due to Aldous [3^. 

We consider the actual value of the pressure. We show that the continuous interpolation method of 
Guerra-Toninelli [T6| applies in a suitably extended form (see [T2] and jH]). We then prove an extended 
variational principle for the considered model and we use it to obtain a rigorous control of free energy 
bounds. In particular we show that an upper bound is obtained by restricting the optimization scheme 
to hierarchical structures (Ruelle probability cascades) which recovers the replica symmetric solution 
heuristically introduced by the physicists. We moreover show that such an upper bound is exact for 
high enough temperatures. Our analysis applies to an arbitrary number g G N of values for the spin 
variables. 

The paper is organized in three parts. In the first part, dealing with existence and simple bounds, 
the antiferromagnetic q-states Potts model on the Erdos-Renyi random graph is defined (Section [3]) 
and the thermodynamic limit of its pressure is proved to exist (Section [3]) . By an estimate using the 
Cauchy-Schwarz inequality, it is also shown in Section |4] that the pressure admits an elementary upper 
bound. 

The part on existence and simple bounds is then followed by the part on the extended variational 
principle. The cavity functional is constructed in Section [5] and it is shown to lead to a variational 
expression for the pressure. A class of optimizers, given by the Derrida-Ruelle random probability 
cascades, is discussed in Section [6] In Section [7] it is shown how this class leads to an upper bound 
for the pressure which coincides with the replica symmetric solution (in particular the trivial replica 
symmetric solution, i.e., the one involving a uniform overlap probability measure, reproduce the simple 
estimate of Section |4]) . This part on the extended variational principle is comparable to the analogous 
results for the Sherrington-Kirkpatrick model ([TU [T]). Note that the SK model is defined on the 
complete graph and the spin glass behavor has its origin in the random couplings of random signs. 
Here frustration is produced by deterministic anti-ferromagnetic couplings, while the randomness is 
given by the underlying spatial structure of the Erdos-Renyi random graph. 

The last part deals with the full control of the high temperature region. In Section [8] it is shown 
that for temperatures high enough the quenched pressure admits the trivial replica symmetric expres- 
sion as a lower bound. This result, combined with the result of Section |4l proves that in the high 
temperature region the pressure is given by the trivial replica symmetric formula. Since the replica 
symmetric solution is known to have a negative entropy at low temperature this proves the existence 
of a temperature below which the solution is not valid anymore. We discuss these implications and 
other points in Section ^ wherein we summarize our results and present some further questions for 
consideration. Interestingly, the lower bound in the high temperature region leads to the study of a 
model defined on the complete graph with a Hamiltonian which is the one of the SK model plus the 
antiferromagnetic Potts-Curie- Weiss model. This last model is studied in Appendix [X] by a suitable 
extension of the techniques developed in |15j . 

2 The model 

We study the antiferromagnetic q-states Potts model on the Erdos-Renyi random graph. We consider 
a set of N vertices, each vertex has attached a spin variables which can take g € N values: ai G 
{1,2, ... ,q}. The edges of the graph are constructed using a set of N'^ independent and identical 
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random variables {Jij} (with i,j = 1,. ..,N) having a Poisson distribution with 

m,j] = ^ • (1) 

In the following we will denote by E the expectation with respect to the randomness of the graph. 
Because of the double counting (we allow Jij to differ from Jj^i) we obtain an average vertex degree 
c. Strictly speaking, in the construction of the Erdos-Renyi random graph we should use random 
variables distributed as Binomial with paranietcr c/{2N). However the two choices are equivalent in 
the thermodynamic limit TV ^ oo and we find it more convenient to work with the Poisson setting. 
The Hamiltonian of the model is given by 

N 

HN{<y) = ^ Ji,jH<^i,<^j) (2) 

where S{x,y) denotes the Kronecker delta function. Note that since the Jij are non-negative, the 
model is anti-ferromagnetic. For a given inverse temperature /3, we consider the random partition 
function 

Zn{c,I3) = Yl e-^"-^^^ , (3) 

ae{l,...,9}«- 

the expectation of a spin function / : {1, . . . , q}^ i-)- M with respect to the random Boltzmann-Gibbs 

state 

and the quenched expectation 

(/(a)) = E[ujifia))] . (5) 

As usual in disordered systems, it will be convenient to define the quenched state for a function of 
the spins of an arbitrary number n G N of real copies. This is obtained by considering the product 
Boltmann-Gibbs state and then averaging over the disorder. Namely, for a function f of n spin 
configurations a-(i),...,o-(") we define 

f2(/(aW,...,a("))) = -1- Y /(aW,...,a("))e-'5(^-«^'^')+-+^-(-^'"» (6) 

and then 

(/(aW,...,a("))) = ETOaW,...,a(")))] . (7) 
The main thermodynamic quantity we will study is the quenched pressure per particle 

Pn{c,P) = ^E[lnZjv(c,/3)] . (8) 

and its thermodynamic limit 

p{c,/3) = lim Pn{c,P) . (9) 

An important observable that will appear later is the sequence (for n > 1) of arrays qN{ri,r2, ■ ■ ■ , r„) 
with (ri, r2, . . . , r„) € {1, • • • q}", which represents the generalized multi-overlap between n spin con- 
figurations cr*^^\ . . . , a^"'\ and is defined as 

1 ^ 

9iv(ri,r2,...,r„) = _ ^ <5(af \ ri)(5(af \ ra) • • • <5(al"\ r„) • (10) 
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3 Thermodynamic limit of the pressure 

Theorem 3.1 The pressure per particle defined in is a superadditive sequence, i.e. for all TV, A^i, N2 
with Ni + N 

NpNiP) > NipnAP) + N2PnM ■ (11) 
Therefore the infinite volume limit of the pressure per particle defined in (0) exists and is equal to 



pW) = suppjv(/3) 

N 



(12) 



Remark 3.2 The fact that pn{c,I3) has a finite upper bound, uniform in N, can be trivially proved 
by the Jensen inequality N^^^{\ii Zn) < N^^\iiM(Zn), where the right hand side is the pressure of 
the antiferromagnetic Potts-Curie-Weiss model. (See, for example, iSS^).) 

Proof: The proof is obtained by interpolation. For a partition of the systems of size N into two 
subsystems of sizes A^i and N2 and ioi a. t £ [0, 1] we consider the following independent Poisson 
random variables: 



Poisson ( 

\xi-t) 



We define the interpolating Hamiltonian 



N 



Poisson 
Poisson 



2Ni 
2No 



N 



i,j = l »ij=l i,j=Ni+l 

which induces an interpolating random partition function 



ZN{c,P,t) = X! 

<Te{i,. ..,<?}" 



/3ff„(<T,t) 



the expectation with respect to an interpolating random Boltzmann-Gibbs 

1 



ZN{c,l3,t) 



E /(-) 



ae{U....qy 



and an interpolating quenched pressure 



(13) 



(14) 



(15) 



(16) 



PNic,[3,t) = -E[Zw(c,/3,i)] 



(17) 



Since pAr(/3, 1) = pn{P) and pn{P, 0) = pN^iP) + PNiil^) the first statement of the theorem (Eq. (fTTjl ) 
follows from the fundamental theorem of calculus if one can show that the interpolating pressure is 
monotonically non-decreasing in t. The derivative of the interpolating pressure reads 



dpN{c,l3,t) 
dt 



2 



1 ^ 



— y lnc.t(e-^^(--''^^)) - ^ y In 



NNi 



Wf (e 



(18) 
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where the following identity has been used: for a vector X = {Xi, . . . , X,„) of m independent Poisson 
random variables Xi with parameter Xi (t) and a function / : N™ R 



dt 



E[/(X)] = E 



— (/(-'^Ij ■ • ■ , + 1, . . . , ^m) - /(-'^l, • ■ • , Xi, . . . , Xm)) 



Expression (|T8)) can be further simplified using the identity 



and the Taylor expansion 

One obtains 

dpN{c,/3,t) _ 
dt ~ 



(19) 

(20) 
(21) 



-E 



E 



(1 - e-'^)" / 



N 



J2 



i,3 = l 



AT 

— Y^u;-i6{a.,a,))-j;^ ^ c.r(^(^., a,)) 



iViVi 



(22) 



which can be rewritten, using the definition of the sequence of generalized multi-overlap arrays in (jlOp . 



as 



dpN{c,li,t) 
dt 



-0\n 



^ n 

n— 1 



(23) 



ri , TV — 1 



From this expression one sees by inspection that the derivative is non-negative, employing a standard 
convexity argument as in [TB] . Therefore monotonicity of the interpolating pressure is established and 
super-additivity (Eq. pT]) ) follows. The second claim of the theorem, i.e., the existence of the thermo- 
dynamic limit of the pressure and its realization as a supremum (Eq. (|12p ) is a standard consequence 
of super-additivity (see [H]). □ 



Remark 3.3 The same computation goes through also for the ferromagnetic model with the change 
(3 I— >■ — /?. However 1 — < for /? > and therefore the series in 1123]) has alternating signs and 
monotonicity can not be derived anymore by inspection. We believe however that the interpolation 
is monotone also in the ferromagnetic case, though in the opposite direction. This belief is based on 
two facts. Firstly, pressure sub-additivity for the ferro-magnetic model on the Erdds-Renyi random 
graph has been checked numerically for small system sizes f^. This is in agreement with a monotonia 
behavior of the interpolating pressure. Secondly, and more importantly, the numerical checks (f^I) 
for the Ising case (q = 2) show that for < i < 1 the series in \2S\) is dominated by the first term 
and therefore one would be left with the same interpolating pressure of the Curie- Weiss model which 
is know to be sub-additive. This is indeed rigorously shown at zero temperature in '91. Although the 
ferromagnetic model has been fully solved in J 10^ . it would be interesting to extend the monotonicity 
result to all temperature. 
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P™' (c, /3) = ^ In f 1 - ) + In g. (25) 



4 "Trivial" pressure bound 

As a preliminary result, we show that the quenched pressure is bounded from above by the so-called 
"trivial replica symmetric pressure" , which improves the Jensen bound of the previous section. 

Lemma 4.1 For all (3 > and N > 1, 

Pw(c,/3)<p™^(c,/3) (24) 

with 

z \ q 
Proof: Using p^ . we have that 

, 1 ^ 

= lnc.(e-^*(---^)). (26) 

By Jensen's inequality and the concavity of the logarithm, this is bounded from above by 

<^lnE|^l-(l-e->(^i^ f^^5(a.,a,)^j , (27) 



ij — l ri — 1 \ i—1 / ri—1 



where also (|20p was used. Observe that 

q / N \ 2 g 

• (28) 

Since all the (jjv(?'i) are nonnegative and add up to 1, wc can use Holder's inequality to get 

1 = E 9Jv(r-i) < 



ri = l 

and thus 



^ E '^^(^i) ■ V^' (29) 



1 ^ 1 

^E'5(^-'-.)^-- (30) 

Hence, 

^P.(c,/?)<-ln(l — ). (31) 

The bound in the lemma is then obtained by using the fundamental theorem of calculus and observing 
that 

pw(0,/3)=lng. (32) 

□ 

5 Extended Variational Principle 

We begin with a proposition that is useful in the ultimate definition of the extended variational 
principle. 
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Proposition 5.1 Suppose that N,AI G N are chosen and let fiM he any measure on 
Then 

Pn{c,I3) < GN,M{c,(i,HM) , (33) 

where 



Gn,m{c, 13, ^im) = M (c, /3, /^Af) - G'^Nui^-i P-> Mm) 



(2) 



with 



1 



N M 



E 

Te{l,...,q}J>^ cre{l,...,g} 

where the Kij 's are i.i.d. Poisson random variables with parameter c/M , and 

/ M 

I^M{T)exp -(3 X LijS{Ti,T. 



(2) 1 



1° E 

Te{l,...,q}J> 



where the Lij 's are i.i.d. Poisson random variables with parameters cN/(2M'^). 

Proof: For a t E [0, 1] we consider the following independent Poisson random variables: 



such that 



2M2 



for all appropriate indices i,j. We define 

N N M M 



i=l 3 = 1 



and we also define 



Zn.m{,C, (3, flM, t) 



and 



Then 



Zn,m{c, (3, ^M,t) 



E '^A^(^) E 

e{l,...,g}*-f <Te{l,...,<;}" 



ae{l,...,qV 



1 (^91 

— E[lnZ7v,M(c,^,A*A/,0)] = pn{c, 13) + G^^[^j{c, 13, ^im) • 



(34) 



(35) 



(36) 
(37) 

(38) 

(39) 
(40) 
(41) 



since if t = then Hn.m{(^, t, 0) splits into a summand only depending on a and one only depending 
on T. Furthermore, 

^E[lnZNM{c,f3,fiMA)] = G^?M(c,/3,MAf). (42) 
Moreover, as in the proof of Theorem 13. 11 one can show that 

^ ( -^E[ln Zn,m (c, 13, fiM , ty 



1 ^ 



— no UTi 



X5:in..(e-^^(--^-)) + ^ X ln..,(e"^^(--- 



i=l 3 = 1 



i,j=l 



-E 



E 

n=l 



7V2 



W Af 



M 



X ^;'(5(a„a,))-^XX^r('5(^.,r,)) + ^ 5] ^r(<^(n,r,)) 



1=1 j=i 



(43) 
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This can again be rewritten in terms of generalized multi-overlaps as 
(l_e-/3)" « 



-E 



.n=l 



— — ^ — — X] nt(qj,{ri,...,rn) ~2qNiri,...,rn)qMiri,...,rn) +qMiri,...,rn) 
(1 - e-'^)" 

ri,...,rn=l 



.n=l 



n 



' n 



(44) 



which is obviously nonnegative for all t. This proves that 

Pn{c,I3) + g'^^]j^{c,I3,hm) = -^E[lnZjv,M(c,/3,AtM,0)] 



< -E[lnZ;v,M(c,/3,/XM,l)] = G^;^(c,^,mm). (45) 



In other words, 



(46) 
□ 



In order to accommodate a limit where M approaches oo, we mention an equivalent version of the 
function Gn,m{c, (3, (j-m)- 

Lemma 5.2 Let /(I), /(2), ... be i.i.d. random variables uniformly chosen from {1, . . . , A''}, and let 
J(l), J(2), . . . be i.i.d. random variables uniformly chosen from {!,..., M}. Independently of this, 
let K be a Poisson random variable with parameter cN and let L be a Poisson random variable with 
parameter cN/2. Then 



and 



In ^ hm{t) ^ exp I -/3^5((7/(fe),Tj(fc)) I 

re{l,...,g}" ae{l,...,g}JV \ fc=l / 



In ^ MM(r)exp I -^^^(Tj(2fe_i),rj(2fc)) 

Te{l,...,g}" 



Proof: This follows from a property of Poisson random variables known as Poisson thinning (or 
sometimes called Bernoulli thinning). Previously we had i.i.d., Poisson random variables Kij and L^. 
Now we have just two Poisson random variables K and L, but we have i.i.d., uniform random variables 
7(fc), J{k). The Poisson thinning property refers to the fact that the families 

Kij =#{k<K : {I{k), J{k)) = {i,j)} and Lij = #{k < L : {J{2k - 1), J(2/e)) = , 

are distributed identically to K^j and Lij. The reader may also check this as an exercise. □ 

Corollary 5.3 Suppose that hm is a random measure on {1, ... , q}'^ . Then 

Pn{c,P) < E[Gn,m{c,P,Hm)] , 
where the symbol E denotes the expectation with respect to fiM ■ 
Proof: For each random realization of iim, we have 

Pn{c,I3) < GN,M{c,l3,fiM) , 
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almost surely, according to Proposition 15.11 So the corollary follows by elementary properties of the 
expectation. □ 



The following Proposition may be viewed as the M ^ oo limit of Proposition 15.11 



Proposition 5.4 Let (^q)^x random sequence such that each fia is positive, and the series 

converges, almost surely. Also let [ja^k ■ a G N , k G N) be a random array of elements of {1, . . . ,q}. 
We write C for the measure which describes the joint distribution of 

((Aia)aeN, (Ta,fc : aeN, fc e N)) . 

We say that C is "exchangeable" if 

N , fc e N)) = ((Aia)aeN , (T„,^(fc) : a G N , fc e N)) , 



V 



for every non-random permutation t: of N which moves only finitely many k 's. Here — indicates 
equality in distribution. 

For an exchangeable C, we define 



where 



and 



In^Ma X! exp I -^^(5(cr7(fe),Ta,fe) 

a=l cre{l,...,q}" 



fc=l 



(47) 
(48) 

(49) 



In^^^exp ~/3'^S{Ta,2k-l,Ta.2k) 
. a=l \ k=l 

where E is the expectation over C as well as the random variables {I{k))'kLiJ ^ ^''^^ L, where we 
assume that (I{k)),K,L are all independent of one another and of ((/ia), (tq^^)), and /(I), /(2), . . . 
are i.i.d, uniform on {1, ... , N} , K is Poisson with mean cN and L which is Poisson with mean cN /2. 
Then we have 

Pn{c,P) < Gjv(c, /?,£), 

for every exchangeable C. 

Remark 5.5 The condition of exchangeability is not very restrictive. Given any non- exchangeable 
C, we may obtain an exchangeable law by a standard symmetrization procedure. The advantage of 
assuming that C is symmetric is that it allows us to simplify the expression of Gjq . Otherwise it would 
be more complicated. 

Proof: For any fixed M, consider a random realization of the sequence {pLa)a£n and {ja^k '■ ct € 
N , fc e {1, . . . , M}). Define the random measure piM on {1, . . . , q}^^ defined from this as 



A* A/ 



for each r S {1, . . . , q}^^ , where 1{- • • } represents the indicator function of the condition given by 
{•••}. This is the empirical measure, but where we merely truncate the full sequence (tq^i, rQ_2, ■ . • ) 
to the first M components of the spin. Another useful way to state the same thing is to notice that 
for any non-random function / : {1, . . . , q}^^ M, we have 



rG{l,...,g}" 



(50) 



a=l 
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In fact, it is not necessary that / is non-random, merely that it is independent of {fJ.a)aefi a-nd 
{Ta,k '■ a £ N, fc S N). Note that fiM is a random measure, but according to Corollary 15.31 this still 
gives an upper bound. Specifically, 

PN{c,(i) < E[GN,M{c,|3,^lM)] , 

where the expectation is over the law C for {p-a)aefi a-nd (Ta^k '■ a G N, fc G N), from which iim is 
derived as a measurable function. 

According to Lemma 15. 2[ we may write 



G 



(1) 



1 



In 



'^A^ ("^^ exp ^ (5 (cr/(fc) , rj(fc) ) j 

Te{l,...,q}" (Te{l,...,g}" V fc=i / 



where J(l), /(2), . . . are i.i.d. random variables uniformly chosen from {1, . . . , N}, and J{1), J'(2), . . . 
are i.i.d random variables uniformly chosen from {!,..., M}, and independently of this, X is a Poisson 
random variable with parameter cN. Note that, according to (|50p . we may rewrite this as 

In Y Mm(t) Y exp I -/3^(5(cr/(fc),rj(fe)) j 

TG{l,...,q}" o-G{l,...,q}« \ k=l I 

1 °° ( ^ \ 

= ^lnXl'^" Y exp I -/3^(5(cr/(fe),T„,j(fe)) j . 

a = l o-e{l,...,g}" V / 

But, conditional on the event that J(l), . . . , J(K) are all distinct elements of {1, . . . , Af}, we have 
equality in distribution, 

1 °° ( ^ \ 

— In^^a Yl exp ( -^^(5(cr/(fc),r„_j(fc)) j 

a=l cre{l,...,g}« V fc=l / 

^^In^l^" XI exp -^^(5(cr/(fc),ra,fe)j 

a = l CTe{l,. ..,<?}« \ fe=l / 



(51) 



where we have replaced the random indices J(l), . . . , J{K) by the non-random indices 1, . . . , iiT, be- 
cause we assumed that (t^j i, . . . , r^ jv/) are exchangeable, meaning equal in distribution under finite 
permutations. Note that here we use the fact that K and J(l), •/(2), . . . are independent of {\ia)a&i and 
{ja,k ■ a € N , fee N). Moreover, conditional on the value of K, the probability that J(l), . . . , J{K) 
are all distinct is 

M{M - 1)---{M - K + 1) 

If we take a single realization of K for all M's, then we see that this conditional probability converges 
to 1, pointwise, almost surely. So we are justified in making the rearrangement in (j51l) . with high 
probability. Moreover, conditioning on K, we see that the function on the left hand side of ([51]) is 
bounded in the interval [log(g) — (3{K/N),\og{q)]. This is summable against the distribution of K. 
Therefore, by the dominated convergence theorem, we have 



lim E 

A/-S-00 



(2) 

A similar argument holds for the second term G)^ {c, 13, C). 

Theorem 5.6 For any c and (3, we have the "extended variational principle, 

p{c,f3) = lim inf G'Ar(c, /3, £) 

N^oo C 



□ 



(52) 
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Proof: We view the equality as a concatenation of an upper bound and a lower bound. The upper 
bound is obtained by optimizing the upper bound in Proposition 15.41 and taking the limit N — ^ oo. 
The lower bound is proved combining sub-additivity with Fekete's lemma [22], 



p(c, B) = lini lim — E 



In 



Z 



N+M 



{c,/3) 



using a particular choice of such that 

1 



lim — E 



In- 



(c,/3) 



Zm{c,I3) 



Zm{c,P) 

This is obtained as a weak limit of Boltzmann-Gibbs measures. By definition, 



(53) 



(54) 



1 



-E 



In- 



Z 



M+N 



Zm(c, P) 



ll 
N 



In 



Ere{l,...,g}*-f So- 



e{i,-,9} 



N C 



where, with J,j ^ Poiss 



2(M+W) 
M 



and -ftTij ^ Poiss 



STe{i,...,9}« 



M N 



-I3Hm(t) 



N 



HM,N{T,cr) ^ ^ Jij6{Ti,Tj) + ^^K.ijS{Tj,ai) + ^ Jijd{cri,aj) 



and, with Jy ^ Poiss {■^) 



j=i i=i 



M 



^m(t) = ^ Jij5{n,Tj). 

We split the latter Hamiltonian as follows. 

M M 

Hm{t) = ^ Jjj5(Tj,r,) + ^ Uj6{n,Tj), 



(55) 



(56) 



(57) 



where Li 



Poiss 



(^ 2M{M+N) \ ^"-^ interested in the limit AI — > oo, so we will ignore the term 



J2i 7=1 Jij^i'^ii'^j) in Hm,n{t,<7), because this converges to almost surely. Hence, 



1e 

TV 



In- 



(c,/3) 



Zw(c,/3) 



1e 

N 



In 



Ere{i,...,q}« E<Te{i.....g}" W exp (-/SE^ii E^i ^^J'5(TJ, ct^ 
llTe{i,...,qV' MmW exp (-^E»^=i^u^(T-i,Tj) 



where 



M 

f^lii^) = exp I -/3 ^ .hj5{n,Tj 



(58) 



Using Poisson thinning we can rewrite this as 



1e 

N 



In 



Z 



M+N 



Zn{c, P) 



= — E 

N 



In- 



e{i,...,9}^ 



e{i,-,<3 



}« Mm(^) exp (-/3 '^(^"./(/c) , cr/(fc))) 



ErG{l....,g}^-^ A^mW exp (-^Efcfl '5(rj(2fc-l),rj(2fe)) 



where Km - Voiss{MNj^) and La/ - Poiss(M %j^^ff_^^) ), and where /(I), /(2), ... and J(l), J(2), , 
are i.i.d. random variables uniformly chosen on {1, . . . , and {1, . . . , Af} respectively. Since 



M 



hm Km = K ^ Poiss(ciV) 



lim Lm — L ^ Poiss 

Af-i-oo 



we obtain formula (1541) with = lim 



lAf^oo fJ-M- 



□ 
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6 Derrida-Ruelle Random Probability Cascade 



Theorem ()5.6|) shows that the cavity functional in eq. (|47| needs to be optimized over measures £. The 
optimal choice of the measure has been conjectured to be related to Derrida-Ruelle random probability 
cascades (RPC). In the following we will show how one can obtain increasing level of Replica Symmetry 
Breaking by considering RPC with an increasing number of levels. We will concentrate on the simplest 
realization of that which produces the Replica Symmetric solution. 



6.1 One level RPC 

Given m G (0, 1), let A^, be the measure on (0, oo): dAm.{x) = mx~"'-~^ dx. Suppose that {^1,^2, ■ • ■ } 
are a Poisson point process with intensity measure Am- This has the important property that the set 



has the same distribution as the set 



for i.i.d., positive random variables Ai, A2, . . . . Let 



Co 

Then using properties of the logarithm 

00 

E In^laAa 



As a preliminary step, we will calculate the upper bound from Proposition 15.41 for {/ia} = {(.a}, with 
all the r's being i.i.d., uniform on {1, . . . ,q}. 

We start with the formula for Gj^ , which is easier. The important thing is to condition on the value 
of L. This is because L is random, but it is not independent for the different values of a. Conditioning 
on L, the multipliers 

Xa = exp ^-/3^(5(Ta,2fc-l,T-a,2fc)^ 

are conditionally independent. We may calculate 



1^ 



k=l 



q q 



So 



E 



In^laCXp -P S{Ta,2k-l,Ta^2k) 



k=l 



L 



L In 1 



1 - e 



Since L has mean cN/2, this gives 



2m 



(59) 



For G^'' we can rewrite this as 



- — F 



N K{i) 



i=l k=l 



12 



where each K{i) = \{k < K : I{k) = i}| is an independent Poisson random variable with parameter 
c. This holds because of the special property for Poisson random variables, that the sum of two 
independent Poisson random variables is Poisson. 

Once again the random variables K{1), . . . , K{N) are the same for all a's, not independent. There- 
fore, we condition on K{1), . . . , K{N). In this case 

(N K(i) 
i=l k=l 



N q 



K{i) 



Jl ^ exp -/3 ^ (5((Tj, Ta,^,k) 



i=l cri = l 



k=l 



The Ta^i^fe's are all independent. In particular they are independent for different values of i. Therefore, 

JV 



E[X^\K{1),...,K{N)] = []e"i— 

i=l 



E 



Using this, if we let k be a Poisson random variable with mean c, then 

' 9 N 



G^,^^ = — E'^lnE^i-'^- 

TO 



E' 



-^Efc=i'5('^.-rfc) 



(60) 



6.2 Trivial bound 

A special case arises in the limit to t 1 • We get 



G 



directly. For the more complicated term, we get 



Gj^^ = E'^lnE^i'-'^" 



.CT=1 



E' 



'"In^E^i--^" L-'^^S^=i*('"'^'=) 

(7=1 

q K 

" 1° E n 



cr=l k=l 

q 



= E«ln 
= E'^ln 



a=l ^ 



q 



In g + cln 1 



q 



Combining this with (|6ip we obtain 



PJv(c,/3) < - ln(^l- ^ 



q 



1 -e- 



In g . 



(61) 
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which matches the upper bound from equation ()24|) . This is not the output of the fuh rephca symmetric 
ansatz. This is the output of the "trivial replica symmetric" ansatz: the one with overlap equal to 0, 
as we will see. 

For the full replica symmetric ansatz, we must also allow for the possibility of a nontrivial overlap, 
which we parametrize by a number p € [0, 1]. This is what we do next. Then the trivial RS ansatz is 
the one with p = 0. 



6.3 Two level RPC 

In order to derive the full replica symmetric ansatz, we need to take a two-level RPC, and then take a 
limiting case. For Proposition l5.4[ all that is required of the /i^'s is that the index set {a} is countable. 
Therefore, we may consider instead of a single countable index, a pair a — (ai, 0:2). 

We let be a Poisson point process with intensity measure A^i. For each ai, we let {£,02(^1)} 

be a Poisson point process with intensity measure . These are independent of each other for 
different values of ai. They are also independent of the first point process {Cqj}- Then we define 



We also define 



We assume that mi < TO2- The reason for this is that Cqi ^ moment only up to (and not 
including) mi, and J2a ^02 ("^i) ^ moment only up to (and not including) m2- But the distribution 
of 



1^ ^(ai,Q2) I 



is the same as the distribution of 



E 



l/mi 



Since we need the sum to be finite almost surely, in order to define ^q,, this means we must have 

7712 > fTT-l- 

In addition to the choice {fJ.a} — {£a}, we also choose an ansatz for the spins Tq,. Actually, we have 
'Ta,k^s, but we will assume that the distribution is i.i.d., for different values of k. So we just focus on 
a single value of k, and suppress that index from the notation momentarily. Let r^^ be i.i.d., uniform 
on {1, . . . , g}. Let also be i.i.d., uniform on {1, . . . , g}, independent for every different choice of 
a — (01,02), even if the ai's coincide (as long as the a2's are different, and vice-versa). Let us also 
choose i.i.d., Bernoulli random variables Xa, with probability p to be 1, and probability 1 —p to be 0. 
Let 

Ta = XaT^^ + {1 — Xa 

In other words, the probability that is just r^^ is p. 
consisting of all a's with the first part ai. 

To calculate G^^\ the easier of the two parts of the cavity field functional, we condition on {^;!,^} 
and on {t^^}. Since the r^'s and Xa's are i.i.d., we do not condition on them. Then we note that 



This is a common value for the "cluster" 



P(r„.2fc-i = Ta,2k I {ri^,k}) = P{Xi =X2 = l)l{ri^ 



2k-i = rl,, 2k} + - P (Xi =X2 = !)]■- 

q 
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This implies that 



E 



1 — e^™^*^ 



For 

we have that 
where 

So we get in general 



Ao = E 



E 



\ _ p-l3S{Tc2k-l,Tc,,2k) 

-m2/3i5(To,2fe-i,-rQ,2fc) | | 



mi/rn2 



-1 1/mi 



InAo = In 

mi 



1\ (l-p2)(i_e-™=^)^^'"^/™^ 1 



1-- 1 ( 1- 



(1 - e-™^^; 



Because the spin fields are independent for different values of k, the effect of the L is just to multiply 
this final answer. Therefore, taking the expectation of that, and dividing by N gives 



mi 



P H 



(1 _ e-™=^) 



mi/m2' 



To get the replica symmetric ansatz, we use the 2-level RPC and take the limits m2 1 1 and mi 4 0. 
Taking m2 t Ij gives 



InAo = — In 

mi 



1 



1 - 



(l-p')(l-e-'^) 



X nil 



q 



p 



Then taking mi | gives 



Thus the replica symmetric value of the first term is 



(2) _ C 
N — c, 



P + 



(1 - e-^) 



(1 - 



{1-e-^)). (62) 



The more complicated term is G^^\ Conditioning on the spins at the first level {r^}, and all the 
K{i) values, we get (using a notation which is clear from the context) 



N 



X^,k},{rlk} 



where the rff. are all i.i.d., uniform on {1, . . . , 5}, and 



q K(i) 

, (7i = l fe=l 



-/3<5(<Ti,ri,fc) 



= Xi^krlk + (1 - ^i,k)Tlk ■ 

The Xi^s are i.i.d., BernouUi-p random variables. Since the formulas are identically distributed for 
different i's and since there are N such z's (canceling the division by N), wc get the formula 



mi 



^{Xk},{Ti} 



q K 



W=l fe=l 



-/3<5(T,Tfc) 
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where once again k is a Poisson random variable with mean c, and now {t^} and {t|} are all i.i.d., 
uniform random variables on {1, . . . , g}, and {Xk} are all i.i.d., Bernoulli random variables with mean 
p, and Tfc = XkT^ + (1 — Xk)T^ for each k. 

We now take m2 t 1 and mi l 0. Taking m2 t 1 gives 



.o-=l fc=l 



which can be rewritten 



-E" InE^^'-'J 



mi 



9 K 



\a=l k=l 









^ mi- 



But 



V 2 



+ {l-p) 1 



Using this and taking the limit mi J, gives 

q K 



l4 E n + (1 - P) (l - ) 



Let us now rewrite this in a manner which is appropriate for taking derivatives at p = 0. We 
can write 6^'^'^'^°'''^'='' = 1 — (1 — e^^)(5(f7, r^). Since the average value of 6{a,Tl) is l/q, we may also 
incorporate that: 



Then we may rewrite E"^*'"^^ |^e-/3'5(o-,Tfc)j 

1 - e-^ 



q 



(l-e-^)(<5(a,ri)-g-i). 



The formula for G^'' is simpler if we introduce a new variable, x = (1 — e ^)/{q — e ^). Therefore, 
we obtain 



G^' = Ing + cln ( 1 - 



1 - e-^ 



+ 



InE" 



n(l-p:^(9<5(^,ri)-l)) 



.k=l 



(63) 



,_(l_e-/i) 



6.4 Local Replica Symmetric stability analysis ot p = and the critical 
temperature 

Now we want to consider this formula as a function of p perturbatively near 0. This allows us to derive 
some results of Zdeborova and Krz§,kala [28 in the present context. We say that the p — RS ansatz 
is "stable to RS perturbations" if it is a local minimizer of the extended variational principle in the 
set of RS ansatze. 



{'2') 

Starting from the simpler term, (I62p . we rewrite G]y as 



In 1 



l-e-^\ , c(g-l) 
2q 



In 1 



p^l-e-P) 



2q 



In 1 



q-{l-e-f^) 



Using X = {1 — e ^)/{q — e ^) this is simpler 

l-e-P 



+ ^ [(? - 1) ln(l + P"^) + ln(l -{q- l)p'x)\ 



(64) 
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This is an even function of p, so only even powers will appear. Taylor expansion shows that 

1 



Therefore, 



This gives 



G 



(q - 1) ln(l + p^x) + ln(l - {q - l)p^x) = 



4^2 



;q{q - l)p X 



In 1 



dp- 



q 



0{p' 



'6c{q — 



(65) 



Now turning to the more difficult term, let us start with 
Then we have 



Let us write f{a, t) — {qS{a, t) — 1) 



G 



(1) 

N 



In g + c In 1 



lnE° 



l[{l-pxfia,Tl)) 



.k=l 



(66) 



As usual, we may interpret the function 



ll{l-pxf{a,Tl)) 



.fc=i 



as a cumulant generating function. But the random variable is multi-linear in p. Therefore, when 
expanding in p, we have to take account of these terms. Also, notice that E"'[f{(j, t)] = E'^[/(cr, r)] = 
as long as the expectations are with respect to the uniform measure. Because of this, various terms 
vanish either in the expectation over E'^ or in the expectation over E^'^^} ^ 

For instance, using the fact that E'^[/(a, r)] = 0, we see that the first derivative in p equals 
0. Moreover, since each factor is linear in p, in taking multiple derivatives (of a single copy of the 
product) means we cannot repeat the derivative of any factor. So we obtain 



dp^ 



■E" 



l[{l-pxf{a,Tl)) 



.k=l 



^ E-[/(cT,rj)/(a,ri)] 



3^k 



But then taking the expectation over K^'^^^ gives because since j ^ k,we have 
E{-'^'>E-[/(a,Ti)/(a,Ti)] = E^ [e^.^ [/(a, r/)] • E^'^ [/(a, r^)]' 



0. 



Continuing, we may easily see that the third derivative is again since E°'[/((T, r)] = 0. Then, the 
next simplest term arises from 



dp^ 



■E" 



\{{l-pxf{a,Tl)) 



^ E-[/(f7, rj)/(a, ri)/(a, r/)/(a, r^)] 



j,k,i,m—l 



Sx" 
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We can rewrite this by expanding the square of the sum, and using rephcated spin variables for products 
of expectations: 



dp' 



■E" 



l[{l~pxfia,Tl)) 



.fc=i 



= ^ E-[/(a,rj)/(a,r,i)/(a,ri)/(a,ri 



j,fc,£,m— 1 



3x4^ J2 E'^^^'[/(a,r/)/(a,ri)/(a',rl)/(a',ri)] 



1 £,m— 1 



Any distinct terms for j, k,£,m vanish in the expectation over E^'^^ ^ . Therefore all must be paired. That 
means that the first summand vanishes entirely. In the second summand, we require {i, m) = {j, k) or 
{£, m) ~ {k,j). These two possibilities give an extra factor of 2. Hence, we obtain 



A calculation gives 



p=0 



J,k=l 



E-[/(a,ri)/(a,ri)] = 



-1 if rj^ri, 

(q-l) if t;=t1. 

Using the i.i.d., uniform distribution on {t^} gives P{t| — t^} = Therefore, 

E{-^'}E-[/(a,Tj)/(a,Ti)] = 0, 
as we claimed before. But now we also have 



1 ifr/^ri, 



which gives 



q-l. 



Therefore, also using the fact that E^y^{{j, k) E {1, . . . , k}^ : j ^ fc}] equals — 1)] — c^, we 



obtain 



'^'^ ^(1) 



-%^{q-\)x^ 



p=0 



(67) 



Since the RS ansatz is cj^-* — G^-^ this means that the local stability of the trivial RS ansatz is 
determined by the sign of the difference between the terms in (j67p and (j65l) 



dp4 



= -6c^{q - l)x'^ + 6c{q - l)x'^ = 6c{q - 1)2;^ [1 - cx'^] . 



We assume q > 1 otherwise the model is trivial (the 1-state Potts model). So the p = solution is 
locally stable if and only if 



1- cx > for 



1 - e- 



q-{l-e~f>) 



Note that (/? = 0) ^ (x = 0) while (/? — )■ 00) =4> (x — s- jzj)- Therefore, we see that the p = solution 
is stable at all temperatures if c < c,, where 



(9-1) 



(68) 
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For instance, for the Ising model which has g = 2, this merely states that c* = 1, which is the 
percolation threshold. For c > c* there is a giant component, and for c < c* there is not. 
If c > c* then the critical temperature for stability is given by the condition 

1- cx = 44> X = —p -i^ -=Vc- 



Since 



1 



q-{l- e-P) 

this means that the critical value of f3 satisfies 

q 



-I ^ Jc ^ 1 



13 _ 



1 -e-f* 



So 



/3* = - In 1 



.-13 - 



1 - 



(69) 



7 Replica Symmetric Solution 

We would now like to compare to Dembo and Montanari's results [TO], as well as those of [11]. So far 
we have imposed the condition that all the r^'s and t|'s are i.i.d., uniform on {1, . . . , g}, and 



Tfc 



if Xk - 1, 



if X. 



0. 



This is a very simple special case of the replica symmetric ansatz. The more general case yields the 
following result. 

Theorem 7.1 The quenched pressure is bounded from above by the Replica Symmetric pressure. 
Namely, 



p{c,P) < Eln 



q K 



En(i-(i 



.-f3- 



ls=l k=l 



-Eln 
2 



l-(l-e-^)^Pi(s)P2(s) 



(70) 



where k is a Poisson random variables with mean c, the Pk = {Pk{s))l^i are i.i.d. random probability 
vectors, i.e. Pk{s) > for each s and X]s=i^fc('*) ~ ^ satisfying the following equality in 

distribution 



(A(s))Li = 



(71) 



Proof: Suppose that J^, P) is a probability space. Suppose that J^i is a sub-cr-algebra such that 
(^ai)a^=i J^i-measurable, and have the appropriate Poisson point process distribution with inten- 
sity C^^d^, as stated before. Also, suppose that (Cqi, 02)02=1 independent of and independent 
of each other for different ai, and distributed according to the Poisson point process distribution 
with intensity for < mi < TO2 < 1- Finally, suppose that the conditional distribution 

{■'"(ai,Q2).fc}a^=iJ conditioned on Ti are independent: 

00 

IP(T(Qi,l),fc e Ai, ^au2),k e ^2, • • • I -T^l) = Jl TP{T(ai,a2),k G I -^1 ) , 



assuming that only finitely many A„'s are different than {1, . . . , g}. Also, assume that the T(Qi,a2),fc's 
are independent for different values of ai. 



19 



Finally, assume that K and L are independent of all the ^^'s and Ta,kS as before. For notational 
simplicity, suppose that there is a such that K and L are Jb-measurable and such that all the ^^'s 
and Ta^'s are independent of J-q. 

Then the extended variational principle works as before with 



G^l^ = — ElnE 
mi 



1 



E 



^1 



-1 mi/m2 



^0 



where we write Tfc for a generic Ta,fe, and k is a Poisson-c random variable, measurable with respect to 
Tq. Similarly, 



— 



-ElnE 



2m 1 

Taking the limit mi 4 and m2 1 1 gives 



E 



-m2/3(5(Tl,T2) 



mi/m2 



G^^ = ElnE 



.s=l 



^1 



and 



G 



(2) 
N 



-ElnE 
2 



g-/35(ri,r2) I jr^ 



The most general choice is to have J^i -measurable i.i.d., random vectors Pk — {Pk{s))l^i chosen 
according to a distribution with support in the set such that Pk{s) > for each s, and J2l=i Pk{s) = 1, 
a.s. Then each Tfc is Pfc-distributed. In this case, 



Gi^) = |Eln 



l-(l-e-^)^Pi(s)P2(s) 



s=l 



and 



Gj^^ = Eln 



q K 



n (1 - (1 - ^"') 



L.s=l /c=l 

Assuming an underlying distribution dp{P), 



l-{l-e-l')Y,P,{s)P^{s) 



s=l 



dp{Pi)dp{P2), 



and 



G^^ = E-^ / In 



^JJ(l_(l_e-V,(a)) 



fe=i 



Taking a Frechet derivative of p in the direction of 5p gives a condition for criticality/extremality: 



In 



l-(l-e-^)^Pi(s)P2(s) 



s=l 



d6p{Pi)dp{P2) 



In 



<7 K 



Ls=i fe=i 



fe=2 



Here we have used the product rule (Leibniz rule) and the symmetry of the formulas with respect 

to permutations of the indices on the Pfc's. In the first integral, both p{dPi) and p{dP2) must be 
differentiated which gives a factor of 2. In the second one, we get an extra factor of k. Note that 
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E[k/(k)] — cE[/(k + 1)]. Note that because of the differentiation, P\ had distribution Jp, and it is 
only the k — 1 last P^'s that had distribution p. Therefore, rewriting P\ as Pq and shifting, we have 



In 



d5p{Po)dp{Pi) 



hi 



g (1 - (1 - e-^)Po(s)) n (1 - (1 - ^~^)P^i')) 



k=l 



dSpiPo) n ^/^(^fc) 



fc=l 



Note that, in order for p + dp to he a probabihty measure, given that p is already a probabihty 
measure, it is necessary that Sp is a signed measure, with total measure 0. Therefore, if we integrate 
any function against dp{dPo), which is constant with respect to Pq, then the integral is zero. Using 
this, and properties of the logarithm, we see that 



In 



l-(l-e-^)£Po(5)Pi(.) 



s=l 



d5p{Po)dp{Pi) 



In 



y (1 - (1 - e-^)Po(.)) nLx(i-(i-e-^)n-Gs)) 



d5p{Po) n dp{Pk] 



k=l 



We have normalized the multiplier, so that summed over s it just gives 0. Therefore, we may rewrite 
this by distributing the two terms in the factor (1 — (1 — e~^)Po(s)), to get 



In 



l-(l-e-^)£Po(5)Pi(s) 



s=l 



d5p{Po)dp{Pi) 



In 



1 n ^ nLl(l^(l-«-^)P^-W) 



fc=l 



Since 5p is general, this implies that the condition for extremality is 



nLi(i-(i-e-^)^fc(^)) 



□ 



8 High temperature region 

A guess at high temperature would be that uJN{S{(Ji, (Jj)) is about 1/q generically (for i ^ j) since that 
is the exact value at infinite temperature. A calculation, using ([25]) and (PB)) . shows that 

Pw(c,/3)-p™^(c,/3) = E[\n{l - ag{(3)uJN,ct,p {5{a,,<j,) - q-^))] dt , 

where ujn,c,i3 is the Boltzmann-Gibbs state at inverse temperature /3, associated to the Hamiltonian 
Hj^{c, a) (here we write explicitly the dependence on the parameter c) and 



1 - (1 - e-^')g-i 



An important observation is that 



a,(/3)[l - q-'] = 1 - e [0, 1) for g > 1 and /3 > 0. 
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By Taylor's theorem, we may write 

n h h 

ln(l - ax) = —7— + ^k+i (a, x) , 



k 
fe=i 



where Rk+i{a,x) is the remainder term, of order (aa;)'^+^: 



We want to consider this expansion with a = aq{f}) and x = 5{ai,aj) — q ^. In particular, we notice 
that for aa; < /(X < 1, we can bound 



-Rfc+i(a,a;) > 



1 / 

1 I ax 



k + 1 \1- n 



If also k is odd, then we know that Rk+i{a, x) > 0. This leads to the family of bounds: for any c, /3 > 0, 
and any integer £ G {1,2,...}, 

N 21-1 k(f>\ i-l 

> p;v(c,/3)-p™^(c,/3) >-^J2E^ E[a,^,,,,,(5(a„a,)-g-i)] dt 

i,5=l k=l 

N „i 



2Af2 
where 

= <(/?)/ 

We will do a second order approximation, so we take £ = 1, and note that 

Jo (l-a9(/3)[l-9 

In terms of doing the second order approximation, since the mean-field antiferromagnet has the correct 
first order behavior and the SK mean-field spin glass has the correct second order behavior, we introduce 
this as an intermediate step. Let us define 



1 ^ 



2N 



which is the mean-field antiferromagnetic term, and 

AT 

f2N 



1 ^ 



where the Yi/a arc i.i.d., standard, normal random variables. All the random variables are supposed 
to be independent. We will henceforth write the original model with a superscript. 



H^^\c,a-) = HN{c,a) , 



in order to look like the other two terms. 
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Let us define 

Zn = Zn{c,X,k) = exp [~I3hI^^ (c, a) - XH^^^ (a) - ^ H^^^ (a) 

o-G{l,...,g}" 

The important parameter /? will be implicit in this notation, since it will stay constant. The parameters 
c, A and k will vary. We write ujn for the expectation with respect to the probability measure on 
{1, . . . , q}^ with weights 



exp 



Then the easiest calculation is 



E 



N 



N 



1 ^ 

-- E[(Ljv(5(a„a,))] 



2iV2 



By Gaussian integration by parts, we have 



(72) 



(73) 



— In Zm 

N 



1 ^ 



But, of course, S^{ai,aj) = 6{cfi,aj). So this gives 

N 



d_ 
Ok 



E 



1 
TV 



InZ 



N 



N N 

^ Yn^N{5{<T,,<r,))]-^ ^E[^2^(^(a„a,))] 



(74) 



For later reference, we note that this implies 

N 

2iV2 



1 ^ 



W7v(<5(CTj,crj))] = ^-Qy^ 



N 



InZ 



AT 



(75) 



and 



1 ^ 9 
- 5:E[^^(<5(a.,a,))] =-2 7r:IE 



2iV2 



= -2 



9k 
9k 



— In Zm 

N 



1 ^ 

— ^ E[(:;^(,5(a„a,))] 



TV 



■InZ 



N 



27V2 

1' 



N 



InZ 



(76) 



Finally, by the Poisson differentiation formula (cfr (126 



dc 



■E 



■InZ 



N 



1 ^ 

— ^ E[ln(l-[l-e-'3]^^(5(a„a,)))] ■ 



2Ar2 



But the power series and the bounds for the logarithm are still valid: 



d_. 
dc 







— In Zm 







1 - e 



-/3 



27V2 



q 

N 



27V2 



N 



Y ^['^N{S{ai,crj) - q ^) 
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Expanding out the terms just involving q ^, we see that 



d_ 

dc 



E 



N 



■InZ 



N 



- 2 \ q ) 2q 



2(j2 



a,(/3)- 262(9, 



N 



&2(g,/3) 

2iV2 



2N^ 

N 



J2 EpAr(<5(a„a,))] 



(77) 



Now we can use the formulas (1751) and ([7S)) to obtain 



9c 



E 







— In Zm 







2g 2(j2 

[a(/3)-26(/3)q-i + fe(/3)]^E 



■InZ 



Therefore, for a fixed /3, we take a curve in (c. A, k) space: 

c(i)=co-i, X{t) ^ [a{f3) - 2b{P)q-'^ + b{P)]t and K(t) ^ 2b{l3)t . 
Then the differential inequality above and the fundamental theorem of calculus gives the bound 



E 



N 



■InZ 



N 



t=ca 



< 



Co 



In 1 



1 - e 



2q 



2g2 



(This is reminiscent of Gronwall's inequality, in that one integrates a differential inequality to get a 
typical bound.) In other words. 



E 



N 



lnZAr(c,0,0) 



> E 



In Zjv(0, [a(/3) - 26(/3)g-i + 6(/?)]c, 26(/3)c) 

1 - e-P 



■In 1 



a{f3) b{P)\ 
2q 2g2 J 



The previous inequality tell us that to study the high temperature of the antiferromagnetic model 
on the Erdos-Renyi we are led to consider Zn{0, X, k) for sufficiently small A and k. Interestingly, 
this is a model we have not seen studied before: it is a mean-field antiferromagnet in conjunction 
with a mean-field spin glass. The two models are conducive to study together since they both satisfy 
inequalities going in the same direction. However, as the mean-field antiferromagnet never has a phase 
transition, it is not surprising to learn that it does not make any significant change to the behavior of 
the spin glass at high temperature. 

Lemma 8.1 For A and k in some open neighborhood of the origin, we have 



lim — ! 

N^oo N 



lnZjv(0,A,K)l = Ing- 

J 2q 4:q^ 



We prove this result in the appendix. For now, we mention how it results in the high temperature 
region for the diluted antiferromagnet. 
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Corollary 8.2 There exist /3* > such that for all /3 < /3* we have the asymptotic lower bound, 



p{c, 13) = lim E 



lnZAr(c,0,0) 



> In (7 + - In I 1 



1-e- 



Combined with Lemma \4.1\ this implies, for all (3 < (3* , 

p{c,P)^p^"''{c,P). 
Proof: (Proof of the corollary.) We already know 



E 



N 



lnZAr(c,0,0) 



> E 



N 



In Z^(0, [a(/3) - 2b{p)q-^ + 6(/?)]c, 26(/3)c) 



-mii- 



2q 2g2 J 



Setting 

and 

we see that 



A = [a{P)-2b{p)q-' +bmc, 
K = 2b{p)c, 



Z^(0,[a(/3)-26(/3)g-i+6(/3)]c,26(/3)c) = Zn{0,\,k). 
But by the lemma, we know that 



lim — E 

N^oo N 

In this case we see that 

A K{q-l) 



In 2'jv(0, A, k) 



Ing 



2q 4g2 



[a{P) - 2b{(3)q-^ + bmc , 2b{l3)c{q-l) 



2q 4q2 



Simplifying, this is 



thus we see that 



hm E 



N 



2q 

aWc biP)c 
2q 2g2 



lnZAr(0, [a(/3) - 2b{(3)q-' + 6(/?)]c, 26(/3)c) 



4g2 



Inq — c 



aW) HP) 
2q 2(^2 



Putting this together with ([7^ gives the corollary. 



(78) 



(79) 



□ 



9 Summary and Outlook 

We have considered the antiferromagnctic Potts model on an Erdos-Renyi graph. While this model 
has some important distinctions from a true spin glass model, it also has many qualitative similarities. 
It is a disordered model in statistical mechanics, and there is a high degree of frustration. We have 
demonstrated that the interpolation method for mean-field and diluted mean-field spin glasses also 
applies to this model. This gives an extended variational principle, analogous to the known results 
for the Sherrington-Kirkpatrick and Vianna-Bray model. In particular, this allows us to interpret the 
replica symmetric ansatz of theoretical physics as a rigorous bound on the free energy. By adapting the 
interpolation method of Guerra and Toninelli, we were also able to prove rigorously that the replica 
symmetric solution is exact for sufficiently high temperatures. 

We now make some remarks about the context of our results, further implications, and important 
open problems. 
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1. It should be mentioned that the second moment method does not work to control the high 
temperature region, because Zfq does not concentrate around its mean. To show this, note that 
for a Poisson random variable J with mean A, we have 

oo » ^ 

Ee-^-^ = ^-^''^YX ^ e"^^''"' = e-^(i-^"'). (80) 
fe=i 

Using this, a simple computation shows that the second moment of Zj^ equals 

nzl\ = exp I -^(1 - e-^) ^ (4(ri) + q^[r,)) 

I ri = l 

+ E 9^(ri,r,)|, (81) 

Compared to 

nZN?= E cxp|-^(l-e-^)E(4(n) + g;^(ri))|, (82) 

cr,(T'e{l,...,g}" I 1-1=1 J 

this grows exponentially faster. 

2. It is known that for low connectivity there is no phase transition. It has been proved indeed [1] 
that the ground state entropy is rigorously equal to the replica symmetric prediction when the 
connectivity is low enough. Because we prove here that for sufficiently high temperatures the 
pressure is given exactly by the trivial replica symmetric formula, this also proves as a subsidiary 
result that the model must have a phase transition for sufficiently large values of the connectivity 
parameter. This is because the entropy of the trivial replica symmetric ansatz 

becomes negative (and physically unrealizable) at sufficiently low temperatures and large enough 
c. 

3. Additionally, due to the extended variational principle, one may prove upper bounds associated 
to any level of replica symmetry breaking, as in |14[ 112) . In particular it has been observed 
pni [55] that there are various transition points, such as a temperature at which the replica 
symmetric ansatz is unstable to a 1-level RSB solution. In our notation, this means that the 
replica symmetry breaking solution is a better optimizer of the extended variational principle. 
This gives a second proof of a phase transition, since their calculations are entirely rigorous for 
the purposes of showing this instability, and one can use an argument as in '24]. We remark 
that this approach implies a phase transition for a temperature at least as low as that required 
to obtain local instability of the replica symmetric ansatz, although, as claimed in I^Dl [2E] for 
various parameter ranges there is a first order transition which happens prior to that point. 

4. The open neighborhood of the origin for the parameters (A, k) referred to in Lemma |8. II is shown 
in Figure [21 where we have relabelled y = X, z = k. Combining this with Corollarv l8.21 one may 
in principle determine explicit bounds for the high temperature region. But these bounds will 
not be sharp. 

5. As an example, we mention that it is known that for each value of q, there is a critical value of 
the connectivity c* (g) such that if c < (g) then the replica symmetric ansatz is correct all the 
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way down to zero temperature. We cannot re-derive this low connectivity result. Technically, 
the obstruction is the appearance of the constants a(/3) and b{l3) multiplying c in equation (|79p 
for the parameters A and k. While the constant a(/3) remains bounded as /3 — > oo, the constant 
b{f3) diverges, which effectively rescales c to oo. It would be interesting to find a solution to this 
problem. 



A Proof of Lemma 18.11 

In this appendix, we will carry out the high temperature analysis to prove Lemma 18.11 Our approach 
is reminiscent of the one carried out for the Sherrington-Kirkpatrick model in 15 . 

A.l Spin glass monotonicity principle 

Suppose that we have two Hamiltonians: a non-random Hamiltonian Hq : — > M, and a Gaussian 
spin glass Hamiltonian, which we write as 

fc=i 

where each '■ — >■ M is a non-random interaction, and — (Ji, . . . , Jmjv) &re i.i.d., Af(0,l) 
random variables. What is important is that 

E[H{a,J)] = 0, 

for each a E and there is a covariance which we may write explicitly in this case as 

A/jv 

E[H{a,J)H{a',J)] = ^ /ifc(a)/ifc(a') • 

k=l 

For each t > 0, we define an interpolating Hamiltonian: 

HNiJ,<y,t) = H^j^\a) + ViHN{a,J) + ^E[H]^{a,J)] . (83) 
We define the random partition function 

and the "quenched" pressure 

PN{t) = j^E[\nZN{J,t)] . (84) 



We also define the random Boltzmann-Gibbs measure 
Lemma A.l For t>0, 



ZNiJ,t) 



^ = -^^[{HN{J',<j)rj^N,] , (85) 

where J and J' are independent realizations of the i.i.d., Af{0, 1) noise. In particular pN{t) is mono- 
tone non-increasing. 
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Proof: Suppose that we have a general Hamiltonian 



HN{J,<J,t,u) = Hf{a)+tH^j^\a) + ^HN{J,cj), 



where H^j^^ (a) is another non-random Hamiltonian. Then 



^ In V e-"" 



(J,a,t,u) 



-1e 

TV 



J,N,t,u 



(86) 



where 



Similarly, 



d_ 

du 



- In y e-ff"('7.- 
N ^ 



2N^ 



E 



But now we use Gaussian-integration-by parts to determine that 



1 



E 



1 



fc=i 



-HN{J,<y,t,u) 



2]V7^ E E ^^^^^^ 



-Hn(J ,cr' ,t,u) 



1 

2iV 



k—l CT^Hiv 

This can be written more concisely as 



k=l cren 



E 



_9_ 

9u 



■E 



In Y 



(J,G,t,u) 



where we write 



1 

2N' 



■E 



(E [hUJ, <y)] „ - (E [Hn{J. a)HNiJ, a')]) 



J,N,t,u 



Putting this together with ([55)1 and specializing to 
and setting u = t, we have 



dt 



N ^ 



t,t) 



27V" 



■E 



,[i7Ar(J,a)iJjv(J,fT')])j,jv,, 



Finally, we can make the inner expectation and outer expectations independent by using independent 
realizations of the noise ^7. Doing this, and using the fact that conditional on the probability 
measure associated to (/(cr, u')) j-^jv,t,u makes a and a' independent and identically distributed, we 
obtain the desired result. □ 
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A. 2 Upper and lower bounds for the spin glass plus mean-field interactions 

We will need to introduce various new terms to the Hamiltonian. Therefore, we take the opportunity 
to redefine some terms and introduce some new ones. In completing the proof we will relate these back 
to the original definitions from the last section. Here are the new terms in the Hamiltonian. 

• The mean-field Potts-type Hamiltonian, 

1 ^ 



The mean-field Potts spin-glass. This has Hamiltonian 

^ 7 



for random coupling constants J^j, which are i.i.d., standard, normal random variables, meaning 
— 1 



• A spin-spin overlap term. This has configuration space ft^ x . We denote the pair as (cr, t) 
for (7, T e {1, . . . , g}^. The Hamiltonian is 

Hillia,T) = ^ E [^('^.,^.)-'Z"'] [S{n,r,)-q~'] . (89) 

We observe some simple calculations in the following lemma. 
Lemma A. 2 There are the following relations among these Hamiltonians: 

E[i/|«,,(a)ij|«,,(r)] =i?(l,(a,r) and 

N 

We leave the proof of Lemma IA.2I as an exercise for the reader. The most important step is the 
identity 

[S{a,,a,)-q-']' = {l-2q-')S{a„<j,)+q-\ 
We are interested in the following quantities: 

Nz 

H^iy, z, a) = 2/i?r/#^(a) + i/|S,,(a) - — , (90) 



Hj,ll{a,a) = {l-2q-')H^,f%a) + -^ 



ZN{y,z) 



E' 



(91) 



PA^(y,^) = ^E[lnZw(y,z)] . (92) 

Guerra and Toninelli proved that the thermodynamic limit of this pressure exists in [17]. It is re- 
markable that they are able to prove the existence of the thermodynamic limit even with mean-field 
ferromagnetic interactions, not just antiferromagnetic terms. More recently, other researchers have 
considered this model in certain regimes outside the high-temperature region, for the Ising case, (? = 2, 
to try to establish the Ghirlanda-Guerra identities [B]. We will not need such results, since we restrict 
attention to a high temperature region. 

Guerra and Toninelli's general spin glass monotonicity principle implies the following. 
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Corollary A. 3 For any y e M, z > 



1 — 2q~^ 

PN{y,z)<pN[y 7i z,0] . (93) 



More precisely, 



d ( 1-29-1 



N (94) 



N 



47V2 

where {■)y,z denotes the random Boltzmann-Gibbs average, 

(/(-))..= g /(a)^;^(^ and (/(a,a%,= ^^"'"'^ W^) ' 

Proof: The proof is an application of Lemma lA.ll Suppose that y and z are fixed, and, let us 
define the interpolating Hamiltonians for the definition of Let 

H^{J,a) = i/|^„,(a). 

According to Lemma [4.21 this has variance 

ilW(a) = E[i72,(^,a)] = {l-2q-')H^,T{a) 

We also define the non-random Hamiltonian 

<(-)= (2/-^^|^-)<.r('x). 

Then we see that the interpolating Hamiltonian, 

may actually be written in terms of the Hamiltonian in ()90p as 

HN{a,t) = HN{m,Sit),a) for y(t) = y - (z - f) , 5(t) = 

Defining pN{t) relative to the interpolating Hamiltonian in (IM|) . we see that it related to 

PN{t) = PN{y{t),z{t)) . 

Therefore, we apply ((SS]) . applied to t = z. Again using Lemma [A. 2 1 to get the covariance, this gives 
precisely ([Ml)- □ 



The inequality (|93p is a simple application of Lemma lA.ll But that lemma will always only lead 
to upper bounds. In order to get lower bounds, we use the explicit form of the right hand side of (|M| 
and introduce a term to the Hamiltonian corresponding to this. We define 

H^j^\x,y,z,a,T) = -xH!,ll{a,T)+y[H[fP^ {a)+H[f^''^ {t)] + V^ [H^^^^ (95) 
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for x,y ^M. and z > 0. The reason we take —x is that we want the two spin configurations a and r to 
be ferromagnetically coupled. We wih define the general partition function for all three of these terms 
as 

(<T,r)Gf2"xn« 



and the quenched finite-volume approximation to the "pressure," 

1 



(2)/ \ 

Then, by a calculation like ((M)) . we obtain 



2N' 



■E 



(97) 



(98) 



where now we write {■■■)x.y.z for the random Gibbs measure associated to the partition function 
Zj^ {x,y, z), which is consistent with {■)y.z in the special case x = 0. 



Because of this, one may rewrite (|M1) as 



^^PN[y 



1 - 2q- 



-t,z + t 



(99) 



x=0 



Guerra and Toninelli's next key step is to replace the derivative by a finite difference approximation. 
A general property guarantees that pn{x^ y, z) is convex in x. Therefore 



(2) 



dPN 
dx 



ix,y,z) 



< 



PN^{x,y,z) - p^^\0,y,z) 1 (2) ^ 1 / N 

= -PN[^,y,z) pN{y,z} 

XXX 



for any a; > 0. So we may obtain an inequality from (j99p : 



t,z + tj > -PN[y + 



t,z + t j - -Pn I 2;,2/H ■ 



-t,z + t 



This leads to the following: 
Corollary A. 4 For any x > 0, 

dt ' 



PN y 



t,z + t 



> -- 



(2) I 

Pn 2;,?/ + 



1 - 2q- 



-t,z + t 



(100) 



This corollary is a direct consequence of the previous inequality. The key fact is that now, turning 
the inequality around. 



dt 



PN [y + 



1 - 2q- 



-t,z + t 



< 



(2) r , 

Pn a;, y + 



1 - 2q- 



-t,z + t 



and now Guerra and Toninelli's monotonicity principle in the form of Lemma I A . 1 1 mav be applied to 
yield and upper bound for the right hand side. 



Corollary A. 5 For x,y z > 0, and t e [0, z], 



PN'ix,y,z) < p^j^' [x + z,y 



(2) 



1 - 2q-^ 



z,0 , 



(101) 
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This is proved similarly to Corollary IA.3[ except that now we use the pair, (cr, r) e ^2N, as the 
basic spin configuration. 

Proof: Now we replace N by 2N and consider the pair (cr, r) G ^2N- Let 

Using Lemma rA.2| this has variance 

ijW (a, r)=E [Hl^iJ, a, r)] ^ (1 - 2q-') (fff.^ (a) + i?f/;*^(r)) + 2H^lUa^ r) + ^ . 
This time, we define the non-random Hamiltonian 

r) = -(x + r) + (^2; - «/;*^(^) + i?f/;*^(r)) . 

The interpolating Hamiltonian is 

H2N{(T,T,t) Hf^{a,T)+VtH2N{J,(T,T)+^-H^^^{cT,T). 

In terms of ([95]) this is 

H2N{<J,T,t) = HN{x{t),y{t),z{t),a) , 

for 

1 — 2g~^ 

= X + z — i , = y {z — t) , z{t) = t . 

Defining p2N(t) relative to the interpolating Hamiltonian H2n{o', t, t) as in ([5^ . but with N replaced 
by 27V, it is related to as 

Therefore, we apply (j85p . applied to t = and t — z. This gives the inequality that was claimed. This 
time, we do not calculate the remainder, so we do not need to use the general covariance formula. □ 

A. 3 Conclusion of Lemma [8.11 

We may now establish Lemma TS. II by completing the proof By combining (IMl) . (|100p and (llOip . we 
have 

f l-2g-i \ , , 

PN [y z, 0] < PN [y, z) 



< e'^/^VM V —z. - / v?) { x + z - t.v 



e'^'^PNiy ^z,0)-/ p'j^> [x + z-t,y -^z,0]—dt. (102) 



Note that we may rewrite this as 

pn^ y - - — — 0) ^ PN [y, z) 



X 



{2^ 

But p]y (a;,y, 0) represents the pressure from a true mean-field classical spin system. This may be 
calculated by large deviation techniques. Thus both the lower and upper bounds may be calculated 
this way. When they match, we will be guaranteed to be in the high temperature region. 
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Lemma A. 6 For each u,v ffi,. let us define 



C{u, v;p) = ^ {pab - q + ^ I 



,6=1 



.6=1 



E 

6=1 



Pab In Pa6 ■ 



:,6=1 



where we assume p is in A4^, the set of all probability distributions p — {pab)a b=i {^^ ■ ■ ■ j*?}^- -^"'^ 
any x, y € M, we have 

lim Inp^^ (x, y, 0) = - max C{x, xq~^ - y, p) ■ 
Proof: For a pair (cr, r) G fi^' x fi^, let us define p(cr, = {Pab{(^,t))1 G TM^, where 



1 ^ 

a,6(o-,T) = — ^(5(cri,a)(5(T,,&), 



for a,b — 1, . . . ,q. Let us also define 



We note that 



N N 



6=1 



Therefore, we may write 
1 ^ 



N 



- E 



2N 
N 



i.j=i 

9 



S{ai,a)S{a-j,a)S{Ti,b)6{t.i,b) - q ^y^^S{ai,a)S{a-j,a) - q 



:,6=1 



6=1 



;,6=1 



6=1 



Similarly, 



Va=l 6=1 

Therefore, using (j96p . we have the non-random partition function when = 0: 

gAf£:(a:,a:ij"^-y;p(o-,r)) 



where 



9 I q 

£{u,v;p) ^ [Pab-q'^f + l E 



E ^ 



.6=1 



a, 6=1 \ a—1 

But standard large deviation calculations give 

ln#{(a,r) Gl]f : p(a,T)e A} 



E 

6=1 



Y.^pab-q ^) 



lim 

N->oo 



N 



max pa6 In Pa 

a, 6=1 
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Figure 1: By Lemma I A. 71 we can deduce that the region in grey is a subset of the region of pairs (u, v) 
such that C{u,v] p) has the symmetric state p = p* as an optimizer. 



for any set A C J\4^ which is the closure of its interior. (See for instance, [26].) Then Varadhan's 
lemma implies that 



N 



hm — = max \£(x,xq - y] p) ~ >^ Paboipab 



N 



This is the desired result because C{u, v; p) — £{u, v; p) — b=i Pab Inpoft- 



□ 



We define p* e to be the symmetric state: 

P*ab = fo'^ each a,b ^ 1, . . . , g. 

We say that {u,v) is a "trivial point" of C if p* is an optimizer of C{u,v; p). In other words, {u,v) is 
a trivial point of C if 

max£(M, = 2lnq. 

Note that in the case that there are several optimizers, we still say that {x, y) is a trivial point if one 
of the optimizers is p* . 

Lemma A. 7 (i) If {u,v) is a trivial point of C, then so is {u',v') for all pairs such that u' < u and 
v' < V. 

(ii) If{ui, vi) and (u2, W2) cire trivial points, then (u, v) is also a trivial point for any convex combination 

u = tui + (1 — t)u2 and v — tvi + (1 — t)v2 , for t G [0, 1] . 
(Hi) For q >2, let 

m = ^^ln(g-l), 
q-2 

which is the critical temperature for the q-state Potts ferromagnet. Then {I3{q^),0) and (0,/3(g)) are 
both trivial points. 

In Figure [1] we have sketched a region of pairs (w, v) which are trivial, in the sense defined above. 
Proof: Suppose u' < u and v' > v. Then for any p ^ A4^, we have 



C{u,v;p)-C{u\v';p) 



a,b=l 



E 



b - q 



E 

6=1 



b - q 
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Hence, relative to the symmetric measure p* such that p*^ = q ^, we have 



£(u, v\ p) — C{u, v; p*) — C{u ,v ; p) ~ C{u , w ; p*) + 



E 



.b=l 



a, 6=1 
2 



E 

6=1 



E (^«^ " 



The last two terms are nonnegative. Therefore if 

C{u\v',p) > C{u',v';p*), 

for some p € Ai^, then it follows that 

£{u,v;p) > C{u,v;p*) , 

for that same p. Part (i) of the lemma is the statement of the contrapositive of this result. 
To prove part (ii), note that for (u, v) = t{ui, ui) + (1 — t)(u2, W2), we have 

C{u,v;p) = tC{ui,vi; p) + {I - t)C{u2,V2; p) , 

for every p e Mg. Therefore, if C{u, v; p) > £{u, v; p*), then either 

C{ui,vi:p) > C{ui,vi\p*) , 



or 



C{u2,V2;p) > C{U2,V2; p*) • 



or both. So, once again, the result follows by taking the contrapositive statement. 

For part (iii) first note that when v = 0, then C{u, 0; p) is precisely the large-deviation optimization 
problem for the mean-field r-state Potts model with r = q^. This proves that {u = /3((Z^), u = 0) is a 
trivial point. This follows from the rigorous analysis of the mean-field Potts model. See for instance, 
[19] or [g. 

When u — the model is not equivalent to the q'-state Potts model. But we may show that the 
optimizer is still the same. Given any p S A^^, let us define two measures in A^^, which are measures 
on just the set {!,..., q}, as p^-^^ and p^^\ the marginals 



^Pab and = ^Pab 



6=1 



Then we may define a new "product measure," p £ Mq, such that 

5 , _ .(1) J2) 

P06 — Pa Pb ■ 

Note that the term attached to v/2 in the partition function is 

E(p«-'?-Y+EU'^-<z-T. 

a=l 6=1 

But the process of taking marginals and product measures stabilizes in the sense that the marginals of 
p are: p'"'^^ = p*^^-' and p^^-* = p'^'. So, if w = then the only difference between the function evaluated 
at p and p comes from the entropy term, and we obtain: 



C{0,v-p)-C{Q,v;p) 



[Pq6 ln(pQ;,) - Pa6 ln(pa6)] • 



„6=1 
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But a standard calculation with entropy (see for example, chapter 1 of [27]) shows that 

- X! ^P°-^ lll(/5afc) - Pab lll(pa&)] = ^ Pab In 
a,&— 1 a,??— 1 

and this is always nonnegative. So, given any p € A^^, one can always replace p by p, and £(0,w; •) 
will not decrease. 
But also, 

C{0,v-p) = C{v-p^^^)+C{v-p^^)), 

where for any p € A^q, we have 

= ^j^iPa-q-'f-j^P-^^P-)- 

This is the large-deviation optimization problem for the mean-field g-state Potts model, at inverse 
temperature /3. Also note that there is not restriction on p^^^ or p^^^: given any two such measures, 
we may have started out by choosing p to be the their product measure. Therefore, the optimization 
problem for £(0,w; •) has the same solution as twice the optimization problem for the g-state Potts 
model at inverse temperature v. (If u > then saying the Potts model has a negative inverse tem- 
perature just means we consider the antiferromagnet instead of the ferromagnet, which has no phase 
transition because of "convexity.") In particular, this proves that the critical value oi v \s v = f3{q). If 
V > —f3{q) then p* is still an optimizer for £(0, v; •). Again, this follows from the rigorous analysis of 
the mean- field g-state Potts model, as in [19] or [6]. □ 




Proof of Lemma 18. It Combining (|102p . Lemma [A. 61 and Lemma [A. 71 we see that 

/ r e*^"" 
Inq < liminf pAr(j/, 2;) < limsuppjv(y, -z) < e^'^lng— / (Ing) dx — \nq, 

N^co N^oo Jo X 

as long as {u, v) is a trivial point of C for 

^ -1 l-2g-i _i 1 

u = X + z and v — [x + z)q — y -\ ^ z — xq — y + — z . 

Here it is required that x > 0. But the various bounds we established for the partial derivatives of 
apply to show that it is Lipschitz. Therefore, even if (u, v) is only a trivial point of C for 

u — z and v = -z~y, (103) 

it is still the case that pi\f{y,z) \nq as N ^ oo. From Lemma I A. 71 once again, {u,v) is a trivial 
point if it lies in the generalized triangle with vertices (0, /3(q)), (/3(g'^),0) and (0,— oo). But ()103|) is 
a linear mapping. So this means the condition for {z,y) is that it lies in the generalized triangle with 
vertices (0, — /3(q)), (/3(g^), ^/3(q^)) and (0, oo). This is shown in Figure[2] In this region, we know that 
PN{y,z) has the limiting value Ing. Note that this does include a neighborhood of (0,0): < 2: < e 
and — (5 < 2/ < (5 for some e, (5 > 0, as desired. 

Now, combining this with (|94|) . since pniy, z) is constant in this region, we see that 

1 ^ 

Y,^[{5{a,,a,)-q-XA 

converges to zero, as — ^ 00, almost everywhere in this region. This also means that for almost every 
(y, z) in this region, 

1 ^ p 
Jp X! ^^('^"^j))?''^ ^ "^"^ asA^^oo,, 
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Figure 2: The region in grey is a subset of the region of pairs (y, z) (drawn with the coordinate pairs 
reversed as {z,y) on the ordinate and abscissa) such that pjsf(y,z) converges to its trivial value \nq. 



where — > indicates convergence in probability. We can use this to determine the desired formula for 



In Zi\[{0, A, k] 



When A 



we have the trivial formula 
1 



N' 



■E 



In Zjv (0,0,0) 



Inq . 



(Recall that we will assume c = 0, throughout.) That serves as a starting point. We may calculate 



1 

— E 

N 



In Zn{0, a, k) 







■ 1 


= Ing + 






Jo dt 


_7V 



dt. 



for any differentiable path {X{t),k{t)) such that: A(0) = 0, A(l) = A, k(0) = 0, k(1) = n. 

We can relate the partial derivatives of this to the partial derivatives of pN{y, z). Recall that in 
our definition of Hj^{y,z) we did not exactly match the definition needed to obtain Zn{Q,X,k) from 
Lemma l8Tl We subtracted some trivial terms to obtain pN{y, z) that made the analysis proceed more 
easily. Specifically, we subtracted from S{(Ti,aj) in (|87l) and (|88l) . and we subtracted Nz/{2q^) in 
([TO)) . But these are all constant shifts, not depending on the spin configurations. Therefore, none of 
these changes affects the random Boltzmann-Gibbs measures, because they cancel in the numerator 
and denominator in the definition of the Boltzmann-Gibbs weights. In other words, considering the 
definition oi Cjm{- ■ ■) from ([72l) . we do have 

^n{' ■ ■) = {■■■)y.z if c = , A = y and k, = z . 

The partial derivatives of A^^^E[ln Zjv(0, A, k)] are calculated in (|73l) and ((74)) : 



^InZjv 



N 



^ Epw(<5(a„a,))] 



and 



Ok 



E 



N 



N 



N N 



4iV2 



But since we know that 



1 ^ 

-J ^(5(cr.,aj)) 



iV2 



p -1 



as — ^ oo, , 



for almost every (y, z) in the specified region, this means that 



d_ 
dX 



E 



N 



In Zat 



7V-5-00 



1 

2^ 



and — E 

OK 



N 



In Zat 



N^o o 1 1 

4g 4(72 
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for almost every pair (A, k) in the same region. 

Now there seems to be a smaU problem. We have established convergence of the derivatives for 
Lebesque-a.e. pair (A, k) in the desired region. But we have to integrate over a 1-dimensional curve, 
which has Lebesgue measure zero. (In fact, we could have concluded convergence a.e. with respect to 
a 1-dimensional Lebesgue measure, but it would have been for a curve, with a specification for the 
slope between y and z, which may not have matched our desired curve.) However, we also know that 

the derivatives in question are all uniformly bounded. That implies continuity of E lnZAr(0, A, k) 
with respect to A and k. So we can average a bit over these points, to turn the contour integral into 
an area integral. Then we can reduce the averaging window, and use continuity to conclude that we 
obtain the original pressure. Therefore, using this, we do conclude that for (A, n) in this region we 
have 



lim — E 

N^oo N 



In Zj\j{0, A, k) 



, A 

2q 



4q^ 



as claimed. 



□ 
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